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ABSTRACT 

Radio and X-ray emission from brown dwarfs suggest that an ionised gas and a mag¬ 
netic field with a sufficient flux density must be present. We perform a reference study 
for late M-dwarfs, brown dwarfs and giant gas planet to identify which ultra-cool ob¬ 
jects are most susceptible to plasma and magnetic processes. Only thermal ionisation 
is considered. We utilise the Drift-Phoenix model grid where the local atmospheric 
structure is determined by the global parameters Teff, log((i) and [M/H]. 

Our results show that it is not unreasonable to expect Hq or radio emission to 
origin from Brown Dwarf atmospheres as in particular the rarefied upper parts of the 
atmospheres can be magnetically coupled despite having low degrees of thermal gas 
ionisation. Such ultra-cool atmospheres could therefore drive auroral emission with¬ 
out the need for a companion’s wind or an outgassing moon. The minimum threshold 
for the magnetic flux density required for electrons and ions to be magnetised is well 
above typical values of the global magnetic field of a brown dwarf and a giant gas 
planet. Na"*", K"*" and Ca“'" are the dominating electron donors in low-density atmo¬ 
spheres (low log(g), solar metallicity) independent of Te®. Mg+ and Fe+ dominate 
the thermal ionisation in the inner parts of M-dwarf atmospheres. Molecules remain 
unimportant for thermal ionisation. Chemical processes (e.g. cloud formation) affect¬ 
ing the most abundant electron donors, Mg and Fe, will have a direct impact on the 
state of ionisation in ultra-cool atmospheres. 

Key words: brown dwarfs, planets and satellites: atmospheres, stars: atmospheres, 
plasma, radio continuum: planetary systems, radio lines: planetary systems 


1 INTRODUCTION 

Ultra-cool objects like brown dwarfs and giant gas planets 
have masses below the hydrogen-burning limit of ~ 0.08 Mq 
( e.g. Burrows et al. 2001). Brown dwarfs are born like stars 
by gravitational collapse, however, they have not sufficient 
mass to achieve the required core temperature to provide a 
steady rate of nuclear hydrogen fusion. As a consequence, 
gravitational collapse provides the only energy source for 
most of the brown dwarf’s lifetime. Cooling and contracting 
during their entire life, brown dwarfs cannot compensate the 
radiative losses by thermonuclear processes. Brown dwarfs 
evolve from a young objects with an effective temperature 
Teff « 3000 K that is comparable to late M-dwarfs and sur¬ 
face gravities like giant gas planets to high-gravity objects of 
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log((/)=5.0 with an effective temperature lower than Teff « 
500 K. Brown dwarfs are fully convective objects, and they 
can be fast rotators. Brown dwarfs and giant gas planets 
form clouds in their atmospheres which have strong feedback 
onto the atmospheric structure due to element depletion and 
due to a high opacity (Helling & Casewell 2014 and refer¬ 
ences there in). Observations by Biller et al. (2013), Buenzli 
et al. (2014) and Crossfield et al. (2014), for example, sug¬ 
gest that brown dwarf atmospheres show a patchy cloud 
coverage. Transit spectroscopy from extrasolar planets sug¬ 
gest that giant gas planets are covered in hazes and clouds 
too (Pont et al. 2008; Gibson et al. 2012; Sing et al. 2015). 
Model simulations suggest that cloud formation prevails for 
a large range of effective temperatures up to 2800 K, and 
for metallicities as low as [M/H]= —5.0 (Witte, Helling & 
Hauschildt 2009). A significant volume of these clouds is sus¬ 
ceptible to local discharge events (Helling et al. 2013): large- 
scale discharges in the upper cloud regions and corona-like 
small-scale discharges in the inner, denser part of the cloud. 
These local discharge events may generate atmospheric elec¬ 
trical storms. At the same time, storms may ionise the local 
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gas. Luque et al. (2014) modelled the ionosphere of Saturn 
and Jupiter finding that the atmospheric electrical storms 
may produce up to 10® cm“® of free electrons below the 
ionosphere. 

Radio and X-ray emission from ultra-cool dwarfs are 
now well established (e.g. 2MASSJ10475385-I-2124234 by 
Route & Wolszczan 2012; 2MASS J13153094-2649513AB by 
Burgasser et al. 2013). Berger (2002) observed twelve sources 
in radio, X-ray and emission between the spectral types 
M8 and L3.5. Sorahana, Suzuki & Yamamura (2014) present 
mid-IR AKARI observations which suggest the presence 
of chromospheric activity in brown dwarfs. Schmidt et al. 
(2015) derive a rise in magnetic activity in form of Ha emis¬ 
sion from SDSS spectra from 2% for early M-dwarfs (MO) 
to 88% for early L-dwarfs (LO). 39% of the L dwarfs in their 
sample observed multiple times are suggested to be variable. 
Such observations indicate that an appropriately ionised gas 
is present in the atmospheres of such cool objects that al¬ 
lows chromospheric heating. The prevailing question is how 
much of such a cool, cloud-forming atmosphere needs to be 
ionized that a chromosphere could form. 

It is suggested that the magnetic field strength of 
brown dwarfs is as high as 10®G (e.g. Shulyak et al. 2011). 
Lynch, Mutel & Giidel (2015) suggest field strengths of 
~ 2.5 ... 2.5kG for the surface field on TVLM 0513-46 and 
2M 0746-1-20. A correlation between radio activity and ro¬ 
tation has not been settled for brown dwarfs given the wide 
parameter range occupied by these objects (McLean, Berger 
& Reiners 2012; Antonova et al. 2013). Williams et al. (2014) 
demonstrate that ultra-cool objects do not follow the clas¬ 
sical Giidel-Benz relation where the X-ray and the radio lu¬ 
minosity from F —M stars correlate (Giidel & Benz 1993). 
This relation was shown to persist for solar flares and active 
rotating binaries (Benz & Giidel 1994). The deviation of 
ultra-cool stars from the Giidel-Benz relation beyond than 
approximately M5 may suggest a change in the dynamo 
mechanism that produces the magnetic field in such ultra¬ 
cool objects (Cook, Williams & Berger 2014). If the atmo¬ 
spheric gas can couple to the strong magnetic field in brown 
dwarfs, the kinetic energy carried by large-scale convective 
motions may be transported to the top of the atmosphere 
and released in form of flares, quiescent or quasi-quiescent 
emissions. The flares are a sudden release of magnetic en¬ 
ergy from the deeper convective layers. The quiescent emis¬ 
sion is a continuous emission but at lower energy levels than 
flares. The most likely mechanism to produce this emission 
is electron cyclotron maser emission from highly relativistic 
electrons or plasma emission alternatively, incoherent syn¬ 
chrotron or gyrosynchrotron emission from relativistic elec¬ 
trons (Hallinan et al. 2006). Flares and quiescent emission 
were reported for example by Burgasser & Putman (2005) 
for two late-type M-dwarfs, with a magnetic field strength of 
1 kG. Tanaka, Suzuki & Inutsuka (2014) study the mass loss 
of hot Jupiters through MHD waves. They suggest the for¬ 
mation of a chromosphere by Alfven wave heating. Tanaka, 
Suzuki & Inutsuka (2014), however, observe, that the gas 
does not need to be fully ionised, i.e. degree of ionisation 
fc < 1, for the mechanism to work. Mohanty et al. (2002) 
analysed magnetic Reynolds number, Rm oc /e, for a set 
of model atmospheres to show why the chromospheric Ha 
activity should be low in rapid rotating brown dwarfs with 
Teff 1500K and solar metallicity. The works by Schmidt 


et al. (2015) and Sorahana, Suzuki & Yamamura (2014) sug¬ 
gest that L-dwarfs should have a chromosphere (but with a 
reduced filling factor) despite having low magnetic Reynolds 
numbers. 

Our paper presents a theoretical framework using a set 
of fundamental parameters to analyse the ionisation and 
magnetic coupling state of objects with ultra-cool atmo¬ 
spheres. This paper focuses on late M-dwarfs, brown dwarfs 
and giant gas planets spanning Teg = 1000 . .. 3000K. The 
approach is also applicable to, for example, protoplanetary 
disks. Only thermal ionisation is considered for this reference 
study against which the effect of additional processes (e.g. 
dust-dust collisions, cosmic ray ionisation, Alfven ionisation, 
lighting, photoionisation) can be tested in future works. Our 
investigations utilise results from Drift-Phoenix ID model 
atmosphere simulations (Helling et al. 2008; Witte, Helling 
& Hauschildt 2009; Witte et al. 2011). This allows us to per¬ 
form an extensive reference study across the late M-dwarf, 
brown dwarf and planetary regime on the basis of a con¬ 
sistently calculated model atmosphere grid for a large set 
of global parameters (Teff, log((j), [M/H]). The aim of this 
study is to identify ultra-cool objects that are most suscepti¬ 
ble to plasma processes by itself or that lead to instabilities 
that trigger the emergence of a strong plasma. This study 
does not include any multi-dimensional atmospheric flows 
and resulting mulit-D radiative transfer effects. 

Section 2 describes our approach and the use of the 
Drift-Phoenix model atmosphere results. Sections 3 and 5 
introduce our theoretical frame work in form of basic plasma 
and magnetic parameters, respectively. Our study shows 
that ultra-cool atmospheres are composed of an ionised and 
magnetised gas. The local degree of ionisation varies largely 
amongst the objects and throughout the atmospheres. While 
a late M-dwarf has a considerable degree of ionisation 
throughout the whole atmosphere, the degree of thermal ion¬ 
isation for a L-dwarf is rather low but may well be enough 
to seed other ionisation processes like for example due to 
Alfven ionisation (Stark et al. 2013). In Sec. 3.3.2 we demon¬ 
strate that electromagnetic interactions can dominate over 
electron-neutral interactions also in regions of a very low de¬ 
gree of ionisation. In Sec. 3.3.3 we investigate the relevant 
length scales effected by electrostatic interactions in the gas 
phase. Section 4 (with additional material in Appendix A) 
contains our assessment of the local equilibrium chemistry 
abundances with respect to electron donors species across 
the late M-dwarfs, brown dwarfs and giant planets regime. 

In Sec. 5 we demonstrate that 30%-50% of the atmo¬ 
spheric volume can be magnetically coupled in L-dwarfs for 
a magnetic field strength of 10®G. The atmospheric volume 
with a degree of thermal ionisation above a plasma thresh¬ 
old value of > 10“^ is, however, considerably lower. Our 
results show that it is not unreasonable to expect ultra¬ 
cool atmospheres to emit H^ or even in radio wavelength 
as in particular the rarefied upper parts of the atmospheres 
are affected by electromagnetic interactions over many pres¬ 
sure scale heights despite having low degrees of ionisation. 
Section 6 contains the discussion of our results in view of 
previous publications. Section 7 presents our conclusions. 
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Table 1. Drift-Phoenix atmosphere structures used: 
Group 1: giant gas planets and young brown dwarfs. 
Group 2: dependence on log(g). 

Group 3: dependence on metallicity. 



Teft [K] 

logls) 

[M/H] 

Group 1 

1000-3000 

3.0 

0.0 


1000 



Group 2 

2000 

3.0, 4.0, 5.0 

0.0 


2800 




1000 

Groups 2000 3.0 -0.6, -0.3, 0.0, +0.3 

2800 


2 APPROACH 

We aim to assemble a theoretical framework that allows us 
to assess the plasma and magnetic character in atmospheres 
of objects across the stellar-planetary boundary, namely for 
late M-dwarfs, brown dwarfs and planets. Onr approach is 
not limited to these objects as the plasma parameters used 
are fundamental properties of a gas rather than of a par¬ 
ticular object. We utilise the grid of Drift-Phoenix model 
atmosphere structures in order to quantify the plasma and 
magnetic characteristics. Drift-Phoenix is a combination 
of two complementary codes, Drift and Phoenix (Helling 
et al. 2008; Witte, Helling & Hauschildt 2009; Witte et al. 
2011). The Drift code (Woitke & Helling 2004; Helling et al. 
2008) solves a system of equations that describes the station¬ 
ary dust formation process of mineral clouds (seed forma¬ 
tion, growth, evaporation, sedimentation, element depletion) 
and interaction between the dust grains and gas (Woitke & 
Helling 2003, 2004; Helling & Woitke 2006; Helling, Woitke 
& Thi 2008). Phoenix is a hydrostatic radiative transfer 
model atmosphere code (Hauschildt & Baron 1999) that de¬ 
termines the resultant thermodynamic structure of the at¬ 
mosphere (local gas temperature Tgaa [K], local gas pres¬ 
sure Pgas [bar] and local electronic pressure Pe [bar]) from 
fundamental stellar parameters (effective temperature Ten 
[K], surface gravity log((/) [cms“^] and metallicity [M/H]). 
When combined with Drift it provides a self-consistent at¬ 
mospheric model that takes into account cloud formation 
and its impact on the thermodynamic structure and the re¬ 
sulting spectral energy distribution. 

The Drift-Phoenix atmosphere simulations model the 
kinetic formation of mixed mineral cloud particle made 
of Ti02[s], Al 203 [s], Fe[s], SiOaN, MgO[s], MgSiOsH, 
Mg 2 Si 04 [s] which effects 6 elements (O, Mg, Si, Fe, Al, Ti). 
Mg, Si, and Fe are the most abundant elements after O and 
C in a gas with solar element abundances. 

Using Drift-Phoenix atmospheric models with a 
range of effective temperatures (Ten = 1000 K - 3000 
K), surface gravity (log( 5 r) = 3.0, 4.0, 5.0) and metallic¬ 
ity ([M/H]= -0.6, -0.3, 0.0 -1-0.3), we evaluate the degree of 
thermal gas ionisation, the plasma parameter (Sec. 3) and 
magnetic parameters (Sec. 5). Applying a separate chemical 
equilibrium code, we evaluate the gas-phase composition of 
various atmosphere models to assess if the dominating elec¬ 
tron donating species changes or remains the same (Sec. 4.1). 

We apply a chemical equilibrium routine to calculate 




Pgas [bar] 


Figure 1. (Tgas-Pgas) input structure from the Drift-Phoenix 
model atmosphere grid. The hottest models represent late M- 
dwarf atmospheres or atmospheres of young brown dwarfs. The 
coolest models represent atmospheres in the planetary regime. 
Top: Group 1. Middle: Group 2. Bottom: Group 3. 
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the chemical composition in more detailed than available 
from the standard Drift-Phoenix ontput. Our main inter¬ 
est is to evalnate the results regarding the ions in the gas 
phase. A combination of 155 gas-phase molecnles (including 
33 complex hydrocarbon molecules), 16 atoms, and various 
ionic species were used under the assumption of local ther¬ 
modynamic equilibrium (LTE). The Grevesse, Asplund & 
Sauval (2007) solar composition is used for calculating the 
gas-phase chemistry outside the metal depleted cloud lay¬ 
ers. No solid particles were included in the chemical equilib¬ 
rium calculations but their presence influences the gas phase 
by the reduced element abundances due to cloud formation 
and the cloud opacity impact on the radiation field, both 
accounted for in the Drift-phoenix model simulations. 

We utilize Drift-Phoenix model atmosphere (Tgas, 
Pgas, Pe) structures and the dust element depleted abun¬ 
dances as input for our calculations. We group the Drift- 
Phoenix atmosphere structures into three groups for an eas¬ 
ier presentation of our results (Table 1). These groups, de¬ 
fined by a range of global parameters, cover young and old 
giant gas planets, young and old brown dwarfs and late M- 
dwarfs. Figure 1 shows the thermodynamic profiles (Tgas, 
Pgas) for each of these groups. As the effective temperature 
increases, the thermodynamic gas temperature of the atmo¬ 
sphere increases as well for a given log{g) and [M/H]. The 
effect of dust on the local gas temperature in the atmosphere 
appears as a step-like temperature change (backwarming) in 
the atmosphere models. 


3 BASIC PLASMA PARAMETERS 

In the following section, we lay out the theoretical frame¬ 
work that we use to characterise the plasma of an atmo¬ 
spheric environment with respect to its electrostatic and 
magnetostatic behaviour. This is inspired by the wealth of 
radio observations of substellar objects (Berger 2002; Berger 
et al. 2010; Route & Wolszczan 2012; Burgasser et al. 2013; 
Williams et al. 2014). To understand these observations, 
radio wavelength (quiescent emission) and X-ray emission 
(flares), the atmospheric gas must couple with the back¬ 
ground magnetic field. As a consequence, free charged parti¬ 
cles produced in the atmosphere would be accelerated along 
magnetic field lines and released into the upper parts of the 
atmosphere. A magnetic coupling of the local gas would also 
be required for Alfven waves to develop and potentially con¬ 
tributing to an acoustic heating of a chromosphere also on 
such ultra-cool objects (e.g. Testa, Saar & Drake 2015). We 
note that ideal and non-ideal MHD simulations require a 
certain degree of ionisation to allow Alfven wave heating to 
develop as possible mechanisms for chromospheric heating. 

For a plasma to exist, the gas needs to be ionised. The 
degree of ionisation, /e measures the extent to which a gas 
is ionised, and it is defined as 

/e = , (1) 

where Pgas and pe are the gas and electron pressure respec¬ 
tively, both in [bar]. Once we have determined the degree of 
ionisation of the atmospheric gas depending on the global 
parameters (Tefi, log((;), [M/H]) we evaluate the plasma fre¬ 
quency to investigate where in the atmosphere electromag¬ 
netic interactions dominate over kinetic collisions between 


electrons and neutrals, 

U^pe Rne- (2) 

Upe [rads“^] is the plasma frequency (i.e. the frequency at 
which the plasma reacts to an imposed or perturbed elec¬ 
tric fields), iZne [s~^] is the electron-neutral collision fre¬ 
quency. Only if Eq. 2 is fulfilled, we can expect the ionised 
gas to undergo electromagnetic interactions that could lead, 
for example, to discharge processes. A more refined insight 
about electrostatic interactions influencing the atmospheric 
gas can be gained by determining the length scales beyond 
which the Coulomb force of a charge does not any more ef¬ 
fects its surrounding. On length-scale larger than the Debye 
length, a gas will be quasi-neutral and no electrostatic forces 
will effect the gas behaviour. Hence, 


Ad < T (3) 

with Ad the Debye length and L the typical length scale of 
the considered plasma, both in [m]. Ideally, this would be 
associated with the atmospheric volume where the ionised 
gas can couple to an external magnetic field. 

In the following subsections we define each of the plasma 
criteria and evaluate them for our model atmosphere grid. 
All equations and natural constants are given in SI units. 
All results, however, have been converted into cgs unit for 
an easier representation in the astrophysical context. 


3.1 Plasma Frequency: tUpe Une 

In a plasma, if electrons are displaced from their equilibrium 
position (assuming a uniform, stationary ionic background), 
a charge imbalance is imposed on the plasma, creating a lo¬ 
cal, restoring electric field. Consequently, the electrons try 
to re-establish charge neutrality resulting in them oscillat¬ 
ing around their equilibrium position with a particular fre¬ 
quency called plasma frequency. The plasma frequency is 
defined as. 


^pe 


yeorue J 


1/2 


(4) 


with Ue the electron number density [m“^], e the electron 
charge [C], nie the electron mass [kg]. If the plasma fre¬ 
quency, cUpe, is greater than the frequency of collisions be¬ 
tween the electrons and neutral particles, i^ne [s~^] then, 
long-range electromagnetic collective effects dominant over 
short-range binary interactions (see Fig. 3). The collision 
frequency for neutral particles with electrons is given by 
Rne = UgasngasUth.e, where Uth.e is the thermal velocity of 
electrons given by Uth,e=(fcBT's/m.a)^^^ [m s“^], Ugas the am¬ 
bient gas density and Ugas the collision, or scattering, cross 
section of particles. The latter is assumed to be Ugas = Tr-Cgas 
with Fgas = th 2 as the atmospheric gas in late M-dwarfs, 
brown dwarfs and most likely in giant gas planets is com¬ 
posed mostly of molecular hydrogen, H 2 . Therefore, the 
collision cross section is approximated by Ugas ~ UHa ~ 
TT • = 5.81 • ird (ruj = 1-36 • 10“^° m). If the 

charged particles collide frequently with the ambient neu¬ 
tral gas {ojpe/vne < 1), theu motion will be determined by 
nearest neighbour interactions and not by collective, long- 
range electromagnetic interactions. 
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Pgas [bar] 


Figure 2. The degree of thermal ionisation, /e = Pe/{pe + Pgas) 
as a measure of free charged particles for M-dwarf, brown dwarf 
and giant gas planet atmospheres. The M-dwarf atmosphere is 
easily ionised by thermal processes. Atmosphere of brown dwarfs 
can only be thermally ionised in deeper layers. Top: Group 1. 
Middle: Group 2. Bottom: Group 3. 



10-12 10-10 10 '® 10 '® 10 '^ 10 '^ 10 “ 10 ' 
Pgas [bar] 


Figure 3. Ratio of plasma frequency of the electrons and the 
frequency of collisions between neutral particles and electrons. 
Electromagnetic interactions dominate over electron-neutral in¬ 
teractions if (jJpe/h'ne ^ 1- Top: Group 1. Middle: Group 2. 
Bottom: Group 3. 
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3.2 Debye length: Ad <SC L 


The Debye length, Ad, is the spatial length scale beyond 
which a plasma can be considered quasi-neutral (rig ~ ~ 

n-gas). For length scales less than the Debye length, a test 
charge will experience the influence of the charge imbalance 
inside the Debye sphere. The Debye length, Ad [m], resulting 
from the solution of the Poisson equation for a non-zero 
charge density near test charge, is defined as 


rige^ J 


(5) 


with ks = 1.38 X 10 [JK and eg = 8.85 x 10 
[F m“^]. A plasma is quasi-neutral if 


Ad ^ L . 


( 6 ) 


For an ionised gas region to exhibit plasma behaviour, it is 
required that over the length scale of the region, the elec¬ 
tron number density is high enough that L is greater than 
the Debye length. The typical length scale of the plasma, 
L [m], considered in the literature (e.g. Mohanty et al. 
2002; Tanaka, Suzuki & Inutsuka 2014) is the pressure scale 
height which depends on the local gas properties and varies 
with 1/g. Typical values for the pressure scale height are 
10® ... 10® cm for a brown dwarf with log((?) = 5. Tanaka, 
Suzuki & Inutsuka (2014) base their length scale on the defi¬ 
nition of the Alfven speed that is of the order of the velocity 
of sound (their Eq.l2). Also their approach results in a typ¬ 
ical length scale of the order of the pressure scale height. 
Associate with the Debye length is the number of charges 
inside a Debye sphere, Nd {plasma parameter, App. B). If 
Nd S> 1, the ionised gas exhibits plasma behaviour. 


3.3 Plasma parameters across the star-planet 
regime 

In the following, we evaluate the plasma criteria for late 
M-dwarfs, brown dwarfs and giant gas planet atmospheres. 
All results have been calculated considering thermal ioni¬ 
sation only and compose our reference study against which 
the need for additional ionisation processes can be derived. 
First we examine for which global parameters and pgaa val¬ 
ues, the gas is ionised above the threshold value of /e > 10“^ 
(Sec. 3.3.1). In Sec. 3.3.2 we demonstrate that long-range 
electromagnetic collective interactions of many charged par¬ 
ticles can dominate over short-range binary interactions also 
in regions of a very low degree of ionisation. Recent brown 
dwarf atmospheric investigations have focused on the degree 
of ionization to characterize plasma behaviour, in this pa¬ 
per we consider multiple parameters to gain a more detailed 
characterization, (e.g. Osten et al. 2015). In Sec. 3.3.3 we 
demonstrate for which length scale ultra-cool atmospheres 
will be effected by electrostatic processes and that it is not 
unreasonable to expect ultra-cool atmospheres to emit Ha or 
even in radio wavelength as in particular the rarehed upper 
parts of the atmospheres fulfill plasma criteria easily despite 
having low degrees of ionisation. 


3.3.1 Degree of ionisation by thermal processes, fe 

Figure 2 shows the degree of thermal ionisation evaluated 
for the same models represented in Fig. 1 (Table 1). Guided 
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Figure 4. Debye length versus the local gas pressure assuming 
Te = Tgas- For length scales less than the Debye length, a test 
charge will experience the influence of the charge imbalance inside 
the Debye sphere. The horizontal line indicates the typical length 
scale of the plasma L=10^m consider in this study (Helling et al. 
2011). Top: Group 1. Middle: Group 2. Bottom: Group 3. 
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Figure 5. Comparison Debye length of different astrophysical 
environments Groupl (Table 1) is plotted here. The red lines rep¬ 
resent Tgfr = 3000 K- 2200 K covering M-dwarfs and young L- 
dwarfs. Blue lines represent Tgff = 2000K-1000 K covering late 
L-dwarfs and giant gas planets regime. The Debye lengths for 
Martian’s atmosphere are plotted as red diamond (Duru et al. 
2008) and black circle (Trotignon et al. 2001), for plasma com¬ 
posed of electrons, water group ions and protons as violet triangle 
(Yaroshenko et al. 2011) and a pure electron plasma as light blue 
square symbols (Kremer et al. 2006). The black line represents the 
typical atmospheric length scale, L=10® m (Helling et al. 2011) 

by these results we consider /e > 10“^ to be a threshold 
above which the gas is partially ionised and it may ex¬ 
hibit plasma behaviour. The above choice of a threshold 
value is supported by results from laboratory experiments 
and laboratory plasma devices (e.g. Tokamak; Diver 2001; 
Fridman 2008). For a fluorescent tube, the degree of ioni¬ 
sation is /e ~ 10“® according to Inan & Golkowski (2010). 
Christophorou & Olthoff (2004) showed that at low tempera¬ 
ture (Tgas ~ 300 — 600 K) and low density (10^^ —10^® cm~®) 
the gas is weakly ionised with /e ~ 10“® — 10~^. If the 
density of the charged particles increases towards /a —> 1 
the gas will be fully ionised. For example, a fully ionised 
gas is assumed in ideal MHD calculations. This threshold, 
fa > 10“^, allows us to derive the atmospheric volume that 
can be considered as an ionised gas (Sect. 3.4, Fig. 9). De¬ 
riving such atmospheric volume fractions will enable us to 
compare the results from different plasma criteria (Eqs. 2, 3) 
and to demonstrate that a gas does not need to be fully 
ionised in order to exhibit collective plasma effects. 

• Group 1: changing Teff (Fig. 2, top) 

(log(s)=3.0, [M/H]= 0.0) 

A solar-metallicity M-dwarf with TeU = 3000 K achieves 
fc > 10“^ in almost the entire atmosphere. For cooler at¬ 
mospheres with Teff 2800 K, /e > 10“^ is only reached 
for pgaa > 10“^ bar. The atmospheric fraction that reached 
/e > 10“^ increases with increasing Tgff. 

• Group 2: changing log(g) (Fig. 2, middle) 

(Teff = 1000 K, 2000 K, 2800 K , [M/H]= 0.0) 

Values of varying surface gravity log((;)= 3.0, 4.0, 5.0 are 
studied here. Models with Teff = 2800 K reach /e > 10“^ 


for pgaa > 10“"^ bar; models with Teff = 2000 K for pgaa > 
10“^ bar. A small part of the upper atmospheric reaches 
/e > 10“^ for Teff=2800 K and log((?)=4.0, 5.0 due to the in¬ 
creasing contribution of Ca-|- with outwards decreasing pgaa 
(compare Fig. 6). Models with Teff = 1000 K have only 
a small fraction of the atmospheric gas that reaches the 
fe > 10“^ threshold. This occurs in the deepest layers of 
the atmosphere 1 < pgaa < 10^ bar. The atmospheric frac¬ 
tion that reached fe > 10~^ increases with decreasing log(g) 
at high pgas and with increasing log(( 7 ) at low Pgas- 

• Group 3: changing [M/H]= 0.0 (Fig. 2, bottom) 

(Teff = 1000 K, 2000 K, 2800 K, log(5)) 

Values of varying metallicity [M/H]= - 0.6, -0.3, 0.0, -1-0.3 
are analysed. Models with Teff=2800 K and Teff = 2000 K 
satisfy fe > 10~^ for all values of metallicity for pgaa > 
10“® bar. A small part of the upper atmospheric reaches 
fe > 10“^ for Teff=2800 K due to an increasing electron 
donation from Ca-|-. Models with Teff = 1000 K have only a 
small fraction of the atmospheric gas that can be ionised for 
Pgas > 10“^ bar. 

All models of non-irradiated atmospheres show a de¬ 
gree of ionisation which increases from a minimum 
(Pgas ~ 10“® — 10“® bar) with increasing local gas pressure 
values towards the deeper layers of the atmosphere. At¬ 
mosphere models with Teff ^2800 K, log(p)= 3.0, [M/H] = 
0.0 can reach D > 10“^ only for high pgas (inner parts 
of the atmosphere). Only one model atmosphere achieves 
fe > 10~^ throughout nearly the entire atmosphere (Teff = 
3000K, log((?)= 3.0, [M/H]= 0.0). For atmospheres of 
late M-dwarfs, brown dwarfs and giant gas planets for 
Teff = 1000 ... 3000K for varying log(g) and metallicity, we 
observe that: 

- The hottest model has the highest thermal degree of ioni¬ 
sation, fe - 

- The lowest value of surface gravity causes an increase of fe 
at high Pgas, however, the highest value of surface gravity 
causes an increase of fe at low Pgas- Both trends are for a 
given Teff and [M/H]. 

- The highest metallicity values cause an increase of fe at 
high Pgas compared to the lowest metallicity models. For 
Teff=2800 K and Teff = 1000 K the lowest value of the metal¬ 
licity causes an increase of fe at low Pgas- Both trends are 
for a given Teff and log((?). 

Section 4 investigates which atoms or molecules are the 
most important electron donors in these cold atmospheres, 
and hence, responsible for the values of the degree of ther¬ 
mal ionisation, fe . We note that the degree of thermal ion¬ 
isation will be influenced by the formation of clouds if the 
dominating electron donors are amongst the most abundant 
condensing species. Ca does not fall into this category. 


3.3.2 Dominating electromagnetic interaction 

The criterion tUpe S> Vne is used to derive where in an ultra¬ 
cool atmosphere the long-range, electromagnetic, collective 
interactions of many charged particles dominates over short- 
range binary interactions in a ionised gas of a certain degree 
of ionisation. Figure 3 shows the results of this criterion for 
the three groups of model atmosphere structures (Fig. 1). 
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• Group 1: changing Teff (Fig. 3, top) 

(log(fl)=3.0, [M/H]= 0.0) 

As Teff increases, the range of the Pgas where ctjpe S> r'ne 
increase too. For models with Teff ^ 2200 K the entire at¬ 
mosphere satisfies this criterion; for Teff = 2000 K almost 
the entire atmosphere; for models with Teff = 1800 — 1400 K 
in the uppermost and for the innermost parts of the atmo¬ 
sphere and for Teff = 1200 K only for 10“^ < Pgas < 10“^ 
bar. The model with Teff =1000 K is too cool to fulfill this 
criterion. 

• Group 2: changing log(g) (Fig. 3, middle) 

(Teff = 1000 K, 2000 K, 2800 K , [M/H]= 0.0) 

Models with Teff = 2800 K satisfy this criterion through¬ 
out the whole atmosphere except for log=5.0 at the highest 
pressures. For Teff = 2000 K and log(g)=3.0 almost the en¬ 
tire atmosphere fulfills this criterion; Teff = 2000 K and 
log(5r)=4.0 only in the uppermost and for the innermost 
parts of the atmosphere; for Teff = 2000 K and log(g) = 
5.0 only for pgas < 10“® bar. Models with Teff = 1000 K do 
not satisfy the criterion. 

• Group 3: changing [M/H]= 0.0 (Fig. 3, bottom) 

(Teff = 1000 K, 2000 K, 2800 K, log(5)) 

Models with Teff = 2800 K and all value of metallicity and 
Teff = 2000 K with [M/H]= -1-0.3 satisfy this criterion in 
the whole atmosphere. For Teff = 2000 K and [M/H]= 0.0 
almost the entire atmosphere fulfills this criterion. For Teff = 
2000 K and [M/H]= -0.3, -0.6 only in the uppermost and 
for the innermost parts of the atmosphere, 10~^ > Pgas > 
10“^; models with Teff = 1000 K and [M/H]= -|- 0.3 only 
for 10“^ < Pgas < lO'^. Models with Teff = 1000 K and 
[M/H]= 0.0, -0.3, -0.6 and do not fulfill this criterion for 
any atmospheric gas pressure. 

Figure 3 demonstrates that the collective, long-range 
electromagnetic interactions dominate over short-range bi¬ 
nary interactions in atmospheres of low degrees of ionisa¬ 
tions, i.e. for 2800 ^ Teff ^ 2000 K. As Teff and the metal¬ 
licity increase, cupe Vne is easier fulfilled at high Pgas, 
however, as Teff increases and the metallicity decreases, 
ijJpe is easier fulfilled at low pgas for Teff = 2800 K, 

1000 K. This effect is counteracted by an increase in log(g). 
The lowest value of log( 5 r) causes a decrease of Wpe /Vne in the 
uppermost parts of the atmosphere and an increase in the 
innermost parts. Consequently, long-range, electromagnetic, 
collective interactions of many charged particles do not re¬ 
quire a complete ionisation of the atmospheric gas, and a 
moderate gas ionisation is sufficient. 

3.3.3 Electrostatically effected atmospheric length scales 

The Debye length. Ad, is compared to a typical atmospheric 
length scale of the order of the pressure scale height, L=10^ 
m (Helling et al. 2011). Figure 4 shows how the Debye length 
changes depending on the local atmospheric gas pressure, 
and where Ad ^ T is fulfilled. 

• Group 1: changing Tefj (Fig. 3, top) 

(log(s)=3.0, [M/H]= 0.0) 

The criterion Ad T is fulfilled throughout the whole at¬ 
mosphere for Teff ^ 1800 K; only for Teff = 1600 K a small 
atmospheric gas volume for pgas < 10“^*^ bar cannot reach 





Figure 7. Distribution of electron donor’s abundance over atoms, 
molecules and ions for a warm model of Teff = 2000K, log(g)=3.0 
and solar element composition. The upper left corner contains the 
element considered. 
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logPgas [bar] 

Figure 6. The dominating thermal electron donors for a subset of effective temperatures for log(g)=3,0 and solar element abundances; 
Na+, K+ and Ca+ in the low-density atmosphere for all Tefi; Mg'*" and Fe”*" in the inner parts for Ted ^ 2800 K. An off-set of 3 orders 
of magnitudes is applied between the models for all > 1000-ft' to allow a better comparison. 


this criterion. For models with Ted ^ 1400 K this criterion 
is fulfilled for pgas > 10~® bar. As Ted increases, the range 
of Pgas where Ad <C T increases. 

• Group 2: changing log(g) (Fig. 3, middle) 

(Ted = 1000 K, 2000 K, 2800 K , [M/H]= 0.0) 

Models with Ted = 2800 K and 2000 K satisfy this criterion 
for all surface gravity values and throughout the whole at¬ 
mosphere. For Ted = 1000 K and all surface gravity values 
then Ad <C T is fulfilled only for pgas > 10“"^ bar. 

• Group 3: changing [M/H]= 0.0 (Fig. 3, bottom) 

(Ted = 1000 K, 2000 K, 2800 K, log(5)) 

Models with Ted = 2800 K and 2000 K satisfy this criterion 
for all metallicities and Pgas. Models with Ted = 1000 K 
satisfy this criterion in the inner atmospheric regions only 
where Pgas > lO”'^ bar. 

Figure 4 demonstrates that the Debye length is gener¬ 
ally very large in the upper atmospheric regions throughout 
the whole regime of ultra-cool objects, i.e. late M-dwarfs, 
brown dwarfs and giant gas planets. In the upper atmo¬ 
sphere the electron density is low causing an increasing De¬ 
bye length; whereas deeper in the atmosphere the electron 
number density is high and so the Debye length is relatively 
lower. According to Cravens (1997), at the top of the Earth’s 
ionosphere Ad ~ 1cm for Te « lOOOK and ~ 10^^cm“^, 
compared to a vertical extant of « 300km of the ionosphere. 
For the solar wind at « lAU, Ad ~ 700cm {Te « 10®K, 


rie ~ 10^cm“®) . A comparison of different values of the De¬ 
bye sphere for different astrophysical environments is pre¬ 
sented in Table 2 and in Fig. 5. Duru et al. (2008) investi¬ 
gate the electron density in the upper ionosphere of Mars, 
Trotignon et al. (2001) the interaction of the Martian’s at¬ 
mosphere with the solar wind. Both consider the presence 
of the dust in the plasma environment. Yaroshenko et al. 
(2011) model a plasma composed of electrons, water group 
ions and protons with the presence of photoemission due to 
the UV radiation. Kremer et al. (2006) work with a pure 
electron plasma. 

Section B provides supplementary material about Nd, 
the average number of charges in the Debye sphere. The 
values for Nd are ;§> 10® in the rarefied upper part of 
the atmospheres (pgas < 10~^ bar) for all M-dwarf, brown 
dwarf and giant gas planet model atmospheres investigated 
here. Values for the above quoted Debye length for the 
Earth ionosphere and solar wind are 

jv^aoiwmd _ ]^q 10 parameter Nd indicates that the gas 

is dominated by many long-range interactions of charged 
particles, rather than the short-range binary interactions. If 
any charged particle is close enough to others then, it in¬ 
teracts in a collective behaviour, not only with the closest 
one. This collective behaviour distinguishes a plasma from 
a kinetic gas. Collective behaviour occurs if Nd 1. Fig¬ 
ure B1 shows that Nu S> 1 is satisfied throughout the whole 
atmosphere for model atmosphere structures presented in 
Table 1. A substantially higher Nd is required for cooler at- 
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Ionized gas: fg > 10 ^ 



Teff [K] 


Figure 9. The fraction of the thermally ionized atmospheres volume, Vg^as/^atm, for /e > 10 ^ and for M-dwarf, brown dwarf and giant 
gas planet atmospheres. 


Table 2. Debye lengths for different astrophysical environments. 
Figures provides a comparison to results of these papers. 


Object 

Te [K] 

Tie [cm A£)[c7n] 

References 

Martian’s 

Ionosphere 

5000 

ro“^ 1.5x10^ 

Duru et al. (2008) 
(Mar’s atmosphere) 


Martian’s 

Ionosphere 

347 

3131 

12222 

86665 

9.78 

4.12x10“' 

5.82x10“'' 

3.47x10“® 

1.3 

19 

10^ 

1.1x10^ 

Trotignon et al. (2001) 
(Mar’s atmosphere) 


1.16x10'' 

50 

10® 


Saturn 

3x10'' 

30 

2x10® 

Yaroshenko et al. (2011) 

Orbit 

7x10'' 

10 

5.6x10® 

(electron, water group ions 

Insertion 

1.7x10® 

2 

1.9x10® 

and protons plasma ) 


4x10® 

0.1 

5x10® 



Laboratory Kremer et al. (2006) 

Experiments 46418 7.5x10*^ 1.7 (pure electrons plasma) 


mospheres than in atmospheres for higher Teff if thermal 
ionisation is considered only (i.e. Te = rgas,LTB)- 


3.4 Comparing different plasma criteria 

Radio and X-ray observations from low-mass objects suggest 
that their atmospheres contain enough free charges to con¬ 
stitute a magnetized plasma (Hallinan et al. 2008). From our 
evaluation of the thermal degree of ionisation (Sect. 3.3.1), 
we chose a threshold for the degree of ionisation of /e > 10“^ 
above which an atmospheric gas can be sufficiently ionised 
that it may exhibit plasma behaviour. In a plasma, electron- 
electron interactions dominate over collisions between elec¬ 
trons and neutral particles, ujpe ^ (Sect. 3.3.2). Ad¬ 
ditionally, for a plasma to be considered magnetized, the 
magnetic field must be sufficiently strong that it significantly 
influences the electron and ion dynamics, which we address 
in Sect. 5. 

We now cast the results in terms of atmospheric volumes 
to allow a comparison between the results for different pa¬ 
rameters. Figure 9 and 10 summarise our findings in terms 
of the volume fraction, with Fg*)!), the thermally 

ionised volume of the atmosphere and Fatm the total atmo¬ 
spheric volume. V))aa is derived by calculating the fraction 
of the atmospheric volume for which /e > 10“^ (Fig- 9). 
Figure 10 visualises the atmospheric volume fraction where 
cupe S> r'ne is fulfilled. The atmosphere volume that reached 
fe > 10“^ and satisfied cupe i^ne is affected by the global 
parameters as follows: 
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Dominating long-range interactions: cOpg » v^e 
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Figure 10. The fraction of the thermally ionized atmospheres volume, V"g^as/^atm for ujpe ^ t'ne for M-dwarf, brown dwarf and giant 
gas planet atmospheres. 

- If Teff increases, then the thermally ionised atmospheric 
volume fraction increases for a given log( 5 f) and [M/H]. 

- The Fg'^as/V'atm, that reaches /e > 10“^ and » J^ne, 
increases if log(^) decreases for a given Teff and [M/H]. 

- Higher values of the metallicity, [M/H], cause a larger frac¬ 
tion of the atmosphere volume to have a sufficiently ionised 
gas that large-scale electromagnetic interactions dominate 
over electron-neutral collisions for a given Teff and \og{g). 

The late M-dwarfs have the largest atmosphere volume frac¬ 
tion, V'/i/a/Vatm, that reached /e > 10“^ that is for model at¬ 
mosphere structures with Tgff = 2600 — 3000 K, log((7)=3.0, 

[M/H]= 0.0 and Teff = 2800 K , log(ff)=3.0, [M/H]=-h0.3. 

For late M-dwarfs and brown dwarfs, the atmospheric gas 
satisfies /e > 10“^ only for half of their atmosphere. For 
planetary objects the fraction of the volume that reaches 
/e > 10“^ becomes increasingly small except for those that 
have the highest value of metallicity and the lowest value of 
surface gravity. 

Models with 2200 Tefi ^ 3000 A", log((7)=3.0, 

]M/H]= 0.0 have the largest atmosphere volume fraction 
that reached tJpe ^ Vne for a given log(g) and [M/H]. Models 
with Teff =2800 K, log(g)=4.0, [M/H]= 0.0; Teff =2800 K, 
log(3)=3.0, [M/H]= -0.6, -0.3, 0.0, -fO.3 and Tes = 2000 K, 
log(gr)=3.0, [M/H]= -1-0.3 have the largest atmosphere vol¬ 
ume fraction that reached tUpe ^ Vne as well. In young brown 
dwarfs, the atmospheric gas volume that reaches cUpe ':$> r'ne 
is more than 50%. The atmospheric gas volume that reaches 


tUpe S> r'ne for planetary objects is smaller than for the rest of 
the objects, i.e. for Teff = 1000 K, log(g)=3.0, [M/H]= -1-0.3. 
Our results show that V[//‘g/V'atm (t'-'pe ^ r'ne) is larger than 
V"gi‘a/Ktm(/e > 10'^) for 1000 KsCTeff s;3000 K. A general 
observation is that despite a relatively low degree of ionisa¬ 
tion, large-scale electromagnetic interactions can dominate a 
considerably larger atmospheric volume than a /e evaluation 
would suggest for all ultra-cool objects in our sample. 


4 MOST ABUNDANT THERMAL IONS IN 
LATE M-DWARF, BROWN DWARF AND 
GIANT GAS PLANET ATMOSPHERES 

We investigate the atmospheric gas-phase composition re¬ 
garding the most abundant local gas ions to demonstrate 
which are the dominating electron donors across the star- 
planet regime based on our non-irradiated Drift-Phoenix 
model atmosphere grid. This investigation allows us to un¬ 
derstand which gas-phase species are responsible for increas¬ 
ing the number of free electrons in ultra-cool atmospheres 
and consequently, responsible of satisfying the plasma pa¬ 
rameters given in Eq. 1-3. It also allows us to understand 
how the chemical composition of the gas is linked with the 
dominating electron donors in the gas-phase. In the case of 
cloud-forming atmospheres, the abundance of the element 
depleted by cloud formation has an effect on electron donors 



Teff [K] 
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Figure 8. Same like Fig. 7. Both, Mg and Fe are influenced by 
dust formation which strongly decreases both elements, resulting 
into the localize large kink. The upper left corner contains the 
element considered. H-binding species are shown for comparison. 

at the location of the cloud. This section serve also as ref¬ 
erence for future investigations on ionisation processes and 
their effect on the gas composition like in (Rimmer & Helling 
2013). 

4.1 Dominating ions across the late M-dwarf, 
brown dwarf and planetary regime 

Figure 6 demonstrates for a subset of effective temperatures 
Teff= 1000 ... 3000 K (log(g)=3,0 solar element abundances) 
that K^, Na+, Ca^, Mg"*" and Fe^ are the dominating ther¬ 
mal electron donors. These ionic species are the most sig- 
nihcant contributors to the electron number density from 
thermal ionisation. Species that have sufficiently low hrst 
ionization potentials and sufficiently high atmospheric num¬ 


ber densities will contribute most effectively to the thermal 
degree of ionization. Therefore, K^, Na"*" and Ca^ provide 
the majority of thermal electrons. Figure 6 demonstrates 
that K"*" is the dominating thermal electron donor where 
Pgas < 10“^ bar for all atmospheres except in an M-dwarf 
atmosphere of Teff= 3000 K. The second dominating elec¬ 
tron donor is Na^ from Tefr= 2000 K. Na^ dominates for 
Teff= 3000 K. For increasing gas pressure, pgas > 10“^ bar, 
Na+ and K+ are the dominating thermal electron donor 
for Tefi= 2000 K, 1400 K, 1000 K. For example, Mg”*" pro¬ 
vides most of the electrons in the Tefi= 3000 K model for 
Pgas > 10“^ bar. The detailed results for all model groups 
are summarise in Tables A1 and A2 in the Appendix A. 

Figures 7 and 8 show the distribution of the most im¬ 
portant electron donating elements over atoms, molecules 
and ions. For example, Na^ is the dominating Na-species 
at high temperature, but the atomic Na followed by NaH 
and NaOH contain most of the element Na at lower tem¬ 
peratures. Ca, Mg and Fe are involved in the formation 
of many molecules. In the case of Ca, CaCH, CaOH and 
Ca(OH )2 reach the abundances not much lower or even 
higher than the atomic Ca. Our investigations show that 
it is not sufficient to determine the local degree of ionisation 
based on one prescribed electron donor species. Such an ap¬ 
proach has been chosen in various complex simulations like 
MHD simulations for protoplanet ary disks (e.g. Sano et al. 
2000) or atmospheric circulation models (e.g. Perna, Menou 
& Rauscher 2010). 

The exact amount with which K, Na, Ca, Mg and Fe 
contribute to the local degree of thermal ionisation will also 
depend on the amount of each element that is chemically 
locked in cloud particles. Element depletion by cloud forma¬ 
tion is fully taken into account for Mg and Fe (see Sect. 2) 
in the Drift-Phoenix atmosphere simulations. The more 
extended work by Helling, Woitke & Thi (2008) (Fig. 6) sug¬ 
gests that the effect of cloud formation on the Ca"*" abun¬ 
dance is negligible. Morley et al. (2012) (Fig. 3) show that 
Na 2 S[s] is thermally stable for Tgas < HOOK and KCl[s] for 
Tgaa < 900K, hence Na^ and would be effected by cloud 
formation in a similar temperature window. The solar ele¬ 
ment abundances for Na (6.17), K (5.08) and Cl (5.5) are 
lower than for Ca (6.31) and it may therefore be reasonable 
to expect a similarly negligible effect of cloud formation on 
the abundances of Na"*" and K^. 


4.2 Summary on electrostatic parameters 

Our reference study suggests that: 

- M-dwarfs have almost the entire atmospheric gas ionised, 
brown dwarfs present only an ionised gas for pgas > 10“^ 
bar and giant gas planets present only a small fraction of an 
ionised gas for values of Pgas > 10“^ bar. 

- Collective, long-range electromagnetic interactions of 
electrons dominate over short-range binary interactions with 
neutrals particles increase as TeS increases in the atmo¬ 
spheres of M-dwarfs and brown dwarfs. Giant gas planets are 
to cool to fulfill this criterion for large ranges of Pgas; only 
for Teff= 1200 K, log(g)=3,0, [M/H]=0.0 and Teff= 1000 K, 
log(g)=3,0, [M/H]= -1-0.3 a small fraction of the atmospheric 
gas at most deeper parts of their atmosphere. 

- Ad <C T is fulfilled for M-dwarfs throughout their atmo- 
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spheres. For brown dwarfs only for pgas >10 ® bar and for 
giant gas planets for pgas >10“® bar. 

K"*", Na^, Ca^, Mg^ and Fe"*" are the dominating thermal 
electron donors, however, K"*", Na^ and Ca"*" provide the 
majority of electrons for Teff= 1000...3000 K for log(g)=3,0 
and solar element abundances. In particular, the degree of 
ionisation is low in the upper atmospheres where the abun¬ 
dances of those ions are low. As their abundances increase, 
/e increases as well. Long-range electromagnetic interactions 
dominating over collisions (Eq. 2) and a zero electrostatic 
forces inside the plasma (Eq. 3) reqnire a snfBcient number 
of free charged particles. Any process (cloud formation, CR 
impact) that impacts the element abundance of the dom¬ 
inating electron donors will affect to the electric state of 
the atmosphere, and the potential coupling to a large-scale 
magnetic field. 


5 MAGNETIZED PLASMA PARAMETERS 

ACROSS THE STAR-PLANET REGIME 

In the previous sections we discnssed that for a plasma to 
exists, the gas needs to be ionised to a certain degree (fe > 
10“^). The plasma frequency was used to investigate where 
in the atmosphere, electromagnetic interactions dominate 
over kinetic collisions between electrons and nentrals. 

The Debye length provides insight about the length 
scales on which electrostatic interactions influencing the 
atmospheric gas. We demonstrated that the atmospheric 
volume fraction, suspected to show plasma be¬ 

haviour, varies largely through the M-dwarf to planetary 
regime. 

A plasma is considered magnetised when the motion 
and dynamics of the charged particles are influenced by an 
ambient magnetic field. This requires that the magnetic field 
is of sufficient magnitude that the charged particles can on 
average participate in at least one Larmor orbit before col¬ 
liding with a neutral atom or dust particle. Otherwise, fre¬ 
quent collisions with the ambient neutrals will dominate 
the dynamical evolution of the plasma particles and the 
influence of the magnetic field will be negligible. In some 
instances, because of the differing mass between electrons 
and ions, the electrons can be magnetised while the ions 
are not. For magneto-fluid descriptions of plasmas (such as 
magnetohydrodynamics) both the electrons and ions need 
to be magnetised. Radio flares (Route & Wolszczan 2012), 
X-ray flares (Berger et al. 2010) and quiescent radio emis¬ 
sion (Williams, Berger & Zauderer 2013) have been observed 
in brown dwarfs, Schmidt et al. (2015) conclude that 45% 
of their active L-dwarfs are also variable. This fraction of 
L-dwarfs, which is lower compared to the 60% of active 
M-dwarfs that were found to be variable in their sample, 
lead Schmidt et al. (2015) to speculate about a brown dwarf 
chromosphere as origin for the observed Ha emission and 
variability. These observations suggest that there should be 
a strong magnetic field and considerable coupling between 
the magnetic field and the atmospheric gas, either to directly 
accelerate free electrons or to allow plasma waves to travel 
into the low-density upper atmosphere and deposit their en¬ 
ergy causing a chromosphere to develop. It is interesting to 
note here that old brown dwarfs (log(g)=5.0) have enough 


time available to build up a chromosphere even if the acous¬ 
tic heating rates should be low due to a insufficient magnetic 
coupling of the atmosphere. For young brown dwarfs, rapid 
rotation may favour chromospheric heating by a better cou¬ 
pling to a stronger magnetic field. 

5.1 Cyclotron frequency versus collisional 
frequency cUc » VcoW 

The cyclotron frequency is the angular velocity with which 
charged particles gyrate around the magnetic field line 
(Boyd & Sanderson 2003), cJc.s = Wx.s/rL.s = qsB/nis in 
[rads“®], where nia [kg], qs [C], B [T] and ux,s [nis“®] are 
the mass of species s, with charge q^ and speed perpendic¬ 
ular to the magnetic field ux,s respectively; and |[^| = B is 
the magnitude of the external magnetic flux density present 
in the medium. For a charged particle’s motion to be dic¬ 
tated by a magnetic field, the particle needs to complete on 
average one gyration before a collision with a neutral atom. 
Formally, a magnetised plasma requires 

^c,s r'ns, ( 7 ) 

where v-^a [s~®] is the collision frequency for neutral particles 
with charged species s. From Eq. 7 we obtain the minimum 
value for the external magnetic field flux, |(^| = B, that 
is needed to satisfied this criteria. Applying the definitions 
for uJc,a and v^a, we derive the critical magnetic flux density 
that is required for the dynamics of the charged particle to 
be influenced by the background magnetic field: 


eB 

rria 

> 

^ gas,e^gas Vs 5 


(8) 
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^k^TaV'^ 

(9) 

^gas.e^gas 1 
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y ms ) 


where the collision, or scattering, cross section is Ugas = 
^■T’gas- The atmospheric gas in late M-dwarfs, brown dwarfs 
and most likely also in giant gas planets is composed mostly 
of molecular hydrogen, H 2 . The collision cross section is ap¬ 
proximated as Ugas ~ UHa ~ TT • rH 2 ^ = 5.81 • 10“®° m® 
(rus = 1.36- 10“®° m). 

Taking the electrons and ions as the particles that are 
influenced by an external magnetic flux density, Eq. 9 be- 


comes 




to 
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me 1 
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e ' 

V , 

) 

Bi 

> 

mi / 

^gas^gas 1 

e \ 

'keTi) 

1/2 

1 

, mi J 

1 


Grouping the constants, we rewrite Eqs. 10 and 11 as Be oc 
n.gas(m.eTe)®^®, with Be as the minimum threshold for the 
magnetic flux density to ensure that the electrons are mag¬ 
netised and Bi oc ngas(?u.iTi)®'®®, with Bi as the minimum 
threshold for magnetic flux density required to ensure that 
an ion, i, is magnetised in Fig. 12. The ion masses, mi, are 
taken to be for K^, Na"*" Ca"*", Fe®", and Mg®" according to 
Sect. 4 assuming local thermal equilibrium, Tgas ~Ti «Te. 

Figure 11 shows that for a Pgas < 10° bar and B= 10® 
G and for a Pgas < 10“® bar and B= 10 G , Wce >> Vne 
(horizontal black line) is reached for all model atmosphere 
structures. There is almost no dependence on Teft, log{g) 
and the metallicity. Largest values of cJce/r'ne are reached 
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Pgas [bar] 

Figure 11. Ratio of cyclotron frequency of the electrons and the 
frequency of collisions between neutral particles and electrons. 
Top: Group 1. Middle: Group 2. Bottom: Group 3. 



Figure 12. Magnetic flux density required for electrons, Be 
(lower set of lines), and ions, Bi (upper set of lines), to be mag¬ 
netically coupled to a background magnetic field in the object 
(B= lOG - giant gas planets (GP), B= 10^ G - M-dwarfs (MD), 
brown dwarfs (BD); black horizontal/vertical lines), li B > B-, 
(or B > Be), ajc,s ^'ns is fulfiled and the gas is magnetised by 
the external magnetic field B. 

for B= 10® G representative for M-dwarfs or brown dwarfs. 
Figure 12 shows for which atmospheric gas pressures, Pgas, 
electrons and ions can be magnetised in the atmospheres of 
M-dwarfs, brown dwarfs and giant gas planets. For M-dwarfs 
and brown dwarfs an background magnetic field flux density 
of B= 10® G is large enough to magnetise the charged parti¬ 
cles: for electrons at pgas < 1 bar and for ions at pgas < 10“® 
bar. For giant gas planets, the magnetised part of the at¬ 
mosphere decreases because of a smaller background field 
(B^ 10 G) compared to M-dwarfs and brown dwarfs. For 
electrons this occurs at Pgaa < 10“® bar and for ions at 
Pgas < 10“® bar. Figure 15 summarises the results on mag¬ 
netic coupling in term of the affected atmospheric gas vol¬ 
ume, V^^aa/Fatm, that reach Wc.e > t'ne: 

- If Teff increases, the magnetically coupled volume of an 
atmosphere increases for a given log((/) and [M/Hj. 

- If log(p) decreases, then the magnetically coupled volume 
increases for a given Teft and [M/Hj. 

- Higher values of the metallicity, [M/H], cause an increase 
of Fgaa/Fatm for a givcu Teff and log(p). 

For a fixed value of magnetic flux density, M-dwarfs and 
brown dwarf atmospheres have the largest magnetically cou¬ 
pled volume. Unsurprisingly, a smaller fraction of a giant gas 
planets atmosphere is magnetically coupled when thermal 
ionisation is considered as the only source of gas ionisation. 
However, this fraction can reach 80% also in a planetary 
atmosphere. The fraction of the atmospheric gas volume, 
Fgili/Fatm, that reaches /e > 10“’’ (Fig. 9) and tjpe » r-ne 
(Fig. 10) increases for the same set of global parameters Teff, 
log((;), [M/H] like the atmospheric gas volume that reaches 
LOce ^ r'ne- Figure 15 demonstrates also that a larger at¬ 
mospheric gas volume can be expected to be magnetically 
coupling than the thermal degree of ionisation had initially 
suggested in Fig. 9. This finding is particularly relevant with 
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respect to the effect that the magnetic field geometry might 
have on the detection of the Ha activity signatures and on 
the radio emission: Donati et al. (2008) show that partially- 
convective M dwarfs host non-axisymmetric large-scale mag¬ 
netic fields with a strong toroidal component, while fully 
convective M dwarfs have stronger large-scale field domi¬ 
nated by a mainly axisymmetric poloidal component. A par¬ 
tial ionisation of the magnetically coupled gas does influence 
the magnetic flux density. 


5.2 Magnetic Reynolds Number, R™ 

The previous sections outlined the framework quantifying 
when an ionised gas in a substellar atmosphere behaves like 
a magnetized plasma. The magnetic Reynolds number is an 
easy to utilise measure of a potentially magnetically cou¬ 
pled ionised gas. Within the context of MHD, the magnetic 
Reynolds number (Rm) is the ratio of the convective and 
diffusive terms from the magnetic field induction equation. 
It quantifies whether the MHD plasma is in the ideal or 
resistive regimes. When the magnetic Reynolds number is 
very large (i.e. in the limit of large length scales), the MHD 
plasma is in the ideal MHD regime and the convective term 
has the dominant influence. In this regime the motion of 
the plasma fluid is determined by the magnetic field and 
vice verse. In the resistive MHD regime, the diffusive term 
is important and dissipative processes as Ohmic dissipation 
(Huang & Gumming 2012; Perna, Menou & Rauscher 2010) 
become significant. The magnetic Reynolds number is de¬ 
fined through the induction equation 

^ = V X (m X R) -f (12) 

where |^| = B [T] is the magnetic flux density, u is the flow 
velocity (formed by electrons and ions), and a [Sm“^] is the 
electric conductivity. The magnetic diffusivity, 77 , is linked 
to the conductivity by ry = l/cr in [m^s“^] and represents 
the diffusion of the magnetic field, a measure of the effect 
of collisions between the electrons and the neutral particles 
on the magnetic field. The collisions between the neutral 
particles and the charged particles (electrons or ions) have 
an influence on the diffusivity of the magnetic field. If the 
effect of the collisions is sufficient to displace them away from 
the magnetic lines, the coupling between the magnetic field 
and the fluid may be not effective. Therefore, the diffusion of 
the magnetic field depends on the frequency of the collisions 
between neutral particles and charged particles. 

The magnetic Reynolds number, Rm can be defined as 
the ratio between the relative strength between the diffusive 
term and the advective term of the induction equation. It 
can be used as a measure of the magnetic coupling calling 
the plasma coupled to the magnetic field if Rm ^ 1 , with 


R 


m 


I V X (u X R) I 
77 I V2R I 


Applying a dimensional analyses, Eq. 13 reduces to 


(13) 


vB/L vL 
rjB / ri rj 


(14) 



Figure 13. Decoupled diffusion coefficient, = c^i^en/^pe 
the Ohmic diffusion coefficient, 77ohm = c^^ei/^pe dom¬ 

inating thermal electron donors for Teff = 1000,2800 K, for 
log(g)=3,0 and solar element abundances; K"*", Na"^, Ca"^, Fe"*" 
and Mg'*" (Fig. 6). The Ohmic diffusion coefficient is smaller than 
the decoupled diffusion coefficient. This result suggest that the 
binary interactions between the ions and electrons (r/ohm) 
not significant compared to the binary interactions between elec¬ 
trons and neutral particles (r]^) in the case of thermal ionisation. 
The lines in both temperature sets appear from top to bottom 
in the order ^Ohm,K+J '^Ohm,Na+’ '^Ohm,Ca+^ 

Voh m,F'e+ ' 


where L [cm] is a typical length scale of the plasma over 
which I I = R varies through a hydrodynamic motion of a 
velocity [dtj, can be approximated by L=10® [m] (Helling et 
al. 2011). The diffusion coefficient, 77 , can be approximated 
by 77 « 77d (Fig. 13). ^ 

Therefore, the expression for the magnetic Reynolds 
number is rewritten as 


Rm « 10 ■ Uflow ( --^TTX I • (15) 

The values for the flow velocity chosen are Vfiow= 10“^ 
cms“^ and vflow= 10® cms“^ guided by values of circula¬ 
tion models (Cooper &; Showman 2005; Menou &; Rauscher 
2009; Rauscher & Menou 2013; Heng & Showman 2015; 
Rauscher &; Kempton 2014). Figure 14 represents the mag¬ 
netic Reynolds number for the Drift-Phoenix model atmo¬ 
sphere structures which are comparable to the earlier results 


^ The diffusion coefficient used in Eq. 14 is given by 7] = rj^-bPohm 
being r/d = c^i^ne/^pe the decoupled diffusion coefficient 
and r/ohm = the Ohmic diffusion coefficient and 

^ ^ohm- Both measure the degree of the dominance of the col¬ 
lisions between electrons-neutral particles and electrons-ions re¬ 
spectively over long-range electromagnetic collective interactions. 
It is easy to relate rj^ with Wpe/i^ne (see Fig. 3). If the latter in¬ 
crease, the diffusion coefficient decrease and the magnetic field 
may be generated and transported by fluid motions allowing the 
magnetic energy to be released into upper layers of the atmo¬ 
sphere as radio. X-ray and Ha emissions. 


























16 M. I. Rodriguez-Barrera, Ch. Helling, C. R. Stark and A. M. Rice 


Magnetized electrons: » v^e 



1000 1500 2000 2500 3000 

Teff [K] 

Figure 15. The fraction of the atmospheric volume, Vgag/Vatm, where ojce t-'ne for M-dwarf, brown dwarf and giant gas planet 
atmospheres. M-dwarfs and Brown dwarfs are represented by B= 10® G and giant gas planets for B= 10 G. 


by Mohanty et al. (2002) that were based on Dusty- and 
Cond-Phoenix: 


- Reynolds number increases as Teff increases because of 
the increasing of ne/ugas in globally and locally hotter at¬ 
mospheres (see Sec. 3.3.1). The highest Rm is reached for 
Teff= 3000 K, log(g)= 3.0, [M/H]= 0.0, vflow= 10® cms"^ 

- Lower gravities cause an increased 7?m for a given Teff 
and [M/H] in the inner, high-density part of the atmosphere 
only. However, higher gravities cause an increased Rn, for a 
given Teff and [M/H] in the outer, low-density atmosphere. 

- Higher values of metallically cause a increased of Rm 
at high values of pgas for a given Teff and log(( 7 ). Low- 
metallicity atmospheres with low effective temperatures, i.e 
planet or T- and Y-dwarf atmospheres, have an increasing 
Rm in the outer, low density. This correlates with a drasti¬ 
cally changing gas-phase chemistry as shown in Fig. 6. 


Our results suggest that ideal MHD, where a fully ionised 
gas is assumed, is best suited for models atmospheres with 
Teff ^ 3000 K which includes M-dwarfs and young brown 
dwarfs. For cooler brown dwarfs and planetary regime ob¬ 
jects only a small fraction of their atmosphere can be con¬ 
sidered in ideal MHD. 


6 DISCUSSION 

6.1 Chromospheres on ultra-cool objects 

Observations in radio, soft X-ray and Hq, wavelengths from 
low-mass objects infer that their atmospheres are populated 
with magnetized plasmas. Radio and X-ray emission from 
ultra-cool dwarfs have been well established by different au¬ 
thors (e.g. Berger (2002), Route & Wolszczan 2012, Bur- 
gasser et al. 2013). Sorahana, Suzuki & Yamamura (2014) 
and Schmidt et al. (2015) suggest the presence of chromo¬ 
spheres in brown dwarfs. Sorahana, Suzuki & Yamamura 
(2014) suggest that weakened H 2 O (2.7^m), CH4 (3.3^m) 
and CO (4.6/im) absorption in combination with moderate 
Ha emission could be linked to chromospheric activity. They 
represent a potential chromospheric heating by an increased, 
constant temperature (with pgas = nkT) in the upper atmo¬ 
sphere of their UCM ID model atmosphere which allows a 
considerably better data fit of their observation. Schmidt 
et al. (2015) use a comparable approach by replacing the 
outer atmospheric temperature of BT-settle model atmo¬ 
spheres with a chromospheric temperature structure where 
the start of the chromosphere, a chromospheric break and 
the start of the transition region are used as parameters. 
Metchev et al. (2015) discuss the likely correlation of mag¬ 
netic spots and high-amplitude photometric variability in 
brown dwarfs with low surface gravity values. Williams et al. 
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Figure 14. Magnetic Reynolds number Rm for the three different 
model atmosphere structure groups described in Sec 2. Rm is 
calculated for a flow speed of Uflow = 10® cms“^. If the flow speed 
increases, then Rm increases. Top: Group 1. Middle: Group 2. 
Bottom: Group 3. 


(2014) demonstrate that ultra-cool objects do not follow the 
classical Gtidel-Benz relationship where the radio luminos¬ 
ity increases proportional to the X-ray luminosity in F — M 
stars (Giidel & Benz 1993). The deviation of ultra-cool stars 
in the Giidel-Benz relationship beyond than approximately 
M5 may suggest a change in the dynamo mechanism that 
produces the magnetic held in such ultra-cool objects (Cook, 
Williams & Berger 2014). Another interpretation of radio 
emission is the concurrence of an auroral region (Nichols 
et al. 2012). Speirs et al. (2014) describe a theoretical ap¬ 
proach for cyclotron radio emission from Earth’s auroral re¬ 
gion providing a physical description for the widely used loss 
cone parameterisation (e.g. Osten et al. 2015). Speirs et al. 
(2014) show that the radiation results from a backward-wave 
cyclotron-maser emission process. The radio emission is gen¬ 
erated by electrons following a horseshoe velocity distribu¬ 
tion, instead of a cone, that travel the magnetic held lines 
downward. The backwards traveling waves cause the upward 
refraction of the radiation which will be further enhanced by 
density inhomogeneities. 

If the atmospheric gas is well coupled with the back¬ 
ground magnetic field in brown dwarfs and planets, the ki¬ 
netic energy carried by large-scale convective motions may 
be transported to the top of the atmosphere and released. 
Tanaka, Suzuki & Inutsuka (2014) suggest that energy from 
the convective part of the atmosphere might be transported 
through the upper atmosphere by magneto-convection pro¬ 
cesses and suggest the formation of a chromosphere by 
Alfven wave heating. Mohanty et al. (2002) carried out a 
study of magnetic field diffusivity to explain why the chro¬ 
mospheric Ha activity in brown dwarfs is low in spite of be¬ 
ing rapid rotators. Mohanty et al. (2002) based their work 
on a grid of model atmospheres (mid-M and L dwarfs) in 
a parameter range Teff = 3000 K..1500 K, log(5r)=5.0 and 
[M/H]=0.0. They explained why in these range of Tefi the 
observation of chromospheric levels activity are lower than 
early M-dwarfs, considering mid-M and L dwarfs as rapid 
rotators. In our work, we extend our model atmosphere grid 
until Teff = 1000 K and we include atmosphere structures 
with different values of log( 5 ) and [M/H] (Table 1). A linear 
held diffusion equation, MHD regime and LTE are used in 
both works. Results obtained for Rm as a measure of the 
ideal or resistive MHD atmosphere (see Fig. 14) could be 
incomplete and therefore misleading. According to Eq. 14, 
Rm oc 1 /t;. Our results demonstrate that the regions where 
Rm > 1 satisfy also fe > 10“^. Atmospheres of ultra-cool 
objects could be ionised and treated as an ideal MHD gas 
only at deep layers. Furthermore, tUce/iZne S> 1 measures the 
coupled between the magnetic field and the atmospheric gas 
and it depends, mostly, on the strength of the magnetic field 
(see Fig. 15). Therefore, it is possible to find large volumes of 
the atmospheric gas that are magnetised, but smaller mag¬ 
netised volumes that are strongly ionised. 

Our results further suggest not only that higher ef¬ 
fective temperatures (in agreement with Mohanty et al. 
(2002)) and higher metallicity atmospheres are the best can¬ 
didates for forming a magnetised atmospheric plasma in sup¬ 
port of radio. X-ray and Hq observations in ultra-cool ob¬ 
jects. Also low surface gravity atmospheres fall in this cat¬ 
egory which supports the interpretation by Metchev et al. 
(2015) that high-amplitude photometric variability in L3- 
L5.5 dwarfs can also be related to magnetic spot appearance. 
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While M-dwarfs have been shown to be fully magnetised, 
L-dwarfs and later brown dwarfs have smaller atmospheric 
volume that can be magnetised in an external magnetic 
field. This findings relate to the activity-vs-SpecT results 
in Schmidt et al. (2015) (e.g. their Fig. 6). The threshold 
of Teft = 2300 K given in their Fig.6 is the Mohanty et al. 
(2002) threshold that points out the limit from which models 
Rm > 1 using v= lO"* cms“^ and 10“^ ^ tj ^10^ (convec¬ 
tion zone) being rj the optical length in J band. Our paper 
suggests that this threshold move towards Tefi = 1400 K for 
the same value of flow velocity and same region. This result 
suggests that atmospheres cooler than Teff = 2300 K may be 
susceptible to be magnetized. Another criterion to consider 
a gas magnetised is cUce ^ i^ne (see Fig. 11). For all mod¬ 
els considered in this study (Table 1)) the atmospheric gas 
fulfill cJce ^ Vne for Pgas < 1 bar. Combining both criteria 
Wee ^ Vne and Rm > 1, a large fraction of possible ac¬ 
tive objects are found for Teff=3000 K-1400 K, log(g)=3.0, 
[M/H]=0.0 (Groupl). 

6.2 Ionisation through non-thermal processes 

Our work focuses on determining global parameters to pro¬ 
vide the suitable local atmospheric conditions for a magne¬ 
tised plasma to be present. The results suggest that ultra¬ 
cool atmospheres are susceptible to plasma and magnetic 
processes even if only thermal ionisation processes are con¬ 
sidered and the influence of dust beyond element depletion 
is neglected. An atmospheric plasma regime and magnetised 
gas were found in the where M-dwarf atmosphere, which 
cooler brown dwarfs and planetary objects require high Pgas, 
hence, their magnetised volume is smaller than in M-dwarf 
atmospheres. 

Additional potentially non-thermal ionisation processes 
will enhance the degree of ionisation and increase the local 
volume affected by a magnetic field. Local enhancement can 
result from dust-dust collisions in large cloud areas (Helling, 
Jardine & Mokler 2011) and Alfven ionisation if the local, 
hydrodynamic wind speed is high enough (Stark et al. 2013). 
Electric storms that develop inside an atmosphere effects the 
extent of an ionosphere causing a link between the local ion¬ 
isation processes and the global effects (Luque et al. 2014). 
Irradiation from a host star for close-in exoplanets or in 
white dwarfs - brown dwarf binaries (Casewell et al. 2013) 
will increase the local thermal ionisation globally. Galactic 
cosmic rays increase the number of free charge particles in 
single brown dwarfs. Gosmic rays are effective at ionising 
the upper atmospheric parts, however, the exact amount is 
hard to quantify without extensive chemistry simulations 
(Rimmer & Helling 2013). 


7 CONCLUSIONS 

We present a reference study for late M-dwarfs, brown 
dwarfs and giant gas planet to identify which ultra-cool ob¬ 
jects are most susceptible to atmospheric gas-phase plasma 
processes. Only thermal ionisation is considered for this ref¬ 
erence study and the influence of dust beyond element de¬ 
pletion is neglected. The effect of additional processes like 
cosmic ray ionisation, irradiation, Alfven ionisation, lighting 
can be evaluated against the reference results in this paper. 


Ultra-cool atmospheres with high Tgff, high [M/H] and 
low log((ji) have large fraction of atmospheric volume where 
plasma processes occur, and are therefore the best candi¬ 
dates for radio. X-ray and Hc< emissions. M-dwarfs have a 
considerable degree of ionisation throughout the whole at¬ 
mosphere, the degree of thermal ionisation for a L-dwarf 
is low but high enough to seed other local ionisation pro¬ 
cesses like Alfven ionisation or lightning discharges. Elec¬ 
tromagnetic interaction dominates over electron-neutral in¬ 
teractions also in regions of a very low degree of ionisation in 
most model atmospheres in our sample. The relevant length 
scales effected by electromagnetic interactions in the gas 
phase are larger in low-density regions of any atmosphere. 
The minimum threshold for the magnetic flux density re¬ 
quired for electrons and ions to be magnetised is smaller 
than typical values of the global magnetic field strengths 
of a brown dwarf and a giant gas planet. A considerably 
lower magnetic flux density is required for magnetic cou¬ 
pling of the atmosphere in the rarefied upper atmosphere 
than in the dense inner atmosphere. Na^, K^and Ga"*" are 
the dominating electron donors in low-density atmospheres 
(low log(g), solar metallicity) independent of Teff. Mg^ and 
Fe'*'dominate the thermal ionisation in the inner parts of 
M-dwarf atmospheres. Molecules remain unimportant for 
thermal ionisation. Ghemical processes (e.g. cloud forma¬ 
tion, cosmic ray ionisation) that affect the abundances of 
Na, K, Mg, Ga and Fe will have a direct impact on the state 
of ionisation in ultra-cool atmospheres. 

Our results suggest that it is not unreasonable to ex¬ 
pect ultra-cool atmospheres (M-dwarfs and brown dwarfs) 
to emit Ha or even in radio wavelength as in particular the 
rarefied upper parts of the atmospheres fulfill plasma crite¬ 
ria easily despite having low degrees of ionisation. Our re¬ 
sults therefore suggest that an ionosphere may emerge also 
in brown dwarf and giant gas planet atmospheres, and that 
the built-up of a chromosphere is likely. Both effects will 
contribute to atmospheric weather features and to space 
weather occurrences in extrasolar, planet-like objects. Ultra¬ 
cool atmospheres could also drive auroral emission without 
the need for a companion’s wind or an outgassing moon. 
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APPENDIX A: MOST ABUNDANT THERMAL 
IONS IN M-DWARF, BROWN DWARF AND 
GIANT GAS PLANET ATMOSPHERES 

Varying Tetr and log(p): Table A1 includes model atmo¬ 
spheres where [M/H] is kept at a constant value of 0.0 and Tgfr 
is varied alongside log(p') (see Sec. 4). The ions shown in the ta¬ 
ble are the most abundant; it does not include ions that are very 
prominent but not the most abundant. Note: When there are two 
ions listed (for example, Al-h/Na-h) the second ion occurs at the 
higher pressure. 

Higher values of log(g) correspond to higher overall values of 
pressure, due to the increased surface gravity. 

It is worth noting that there are other ions which are highly 
prominent within these models, however the following tables only 
show which ones are the most prominent. 

Varying Teff and [M/H] Table A2 shows model atmospheres 
where log(p') has been kept at a constant value of 3.0 and Tgff is 
varied alongside [M/H]. 
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Table Al. First, second and third ion represent the most abundant positive ion in the high, middle and lower pressure regions respectively. 
First ion corresponds to the higher pressure (and therefore higher abundances) and visa verse. 


Teff 

log(g) 


3.0 



4.0 


5.0 



1000 

K+ 

K+ 

K+ 

K+ 

K+ 

K+ 

K+ 

K+ 

K+ 


Na+ 



Na+ 



Na+ 



1200 

K+ 

K+ 

K+ 








Na+ 









1400 

K+ 

K+ 

K+ 








Na+ 









1600 

Na+ 

K+ 

K+ 







1800 

Na+ 

K+ 

K+ 







2000 

Na+ 

K+ 

K+ 

Na+ 

K+ 

K+ 

Na+ 

K+ 

K+ 


Ca+ 









2200 

Na+ 

K+ 

K+ 







2400 

Na+ 

K+ 

K+ 








Mg+ 









2600 

Mg+ 

K+ 

K+ 









Na+ 








2800 

Mg+ 

K+ 

K+ 

Na+ 

K+ 

K+ 

Na+ 

K+ 

K+ 



Na+ 


Mg+ 

Na+ 


Mg+ 



3000 

Na+ 

K+ 

Na+ 








Mg+ 

Na+ 

K+ 








H+ 
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APPENDIX B: PLASMA PARAMETER: 
NUMBER OF PARTICLES IN A DEBYE 
SPHERE: No > 1 

A plasma has the capacity to screen a single charged particle 
placed at any point. That means that any single charged particle 
attracts oppositely charged particles producing a screening and 
repels those who have the same charge. A net space is produced 
in the neighbourhood of any single charged particle, reducing the 
electric field generated by it. The effective range of the net force 
between particles is restricted to the order of the Debye length 
(see Sec. 3). As a consequence, a test particle in the Debye sphere 
interacts only with particles that lie within this sphere. A’d mea¬ 
sured the efficiency of this screening and allows us to calculate 
how many gas particles are required to participate. Hence, only 
particle inside this screening areas (Ad) can be considered as elec¬ 
trostatically active. 

All particles have a thermal velocity due to the tempera¬ 
ture of the plasma. The deflected angle due to the electrostatic 
interactions is bigger if the number of particles around of the 
screened particle in the Debye sphere is small. The movement of 
the screened particle will not be smooth, unlike in the situation 
when the number of particles in the Debye sphere is sufficiently 
large to reduce it. That is why the Debye length increases as 
the number of particles in the screened sphere decreases. This 
is demonstrated in Fig. 4 where all Debye length increase with 
height in the atmosphere, i.e. with the outward decreasing local 
gas pressure. 

The change in velocity due to the interactions with the par¬ 
ticles produces a non-negligible net electrostatic force inside the 
Debye sphere. Therefore, large numbers of particles that are uni¬ 
formly distributed inside the Debye sphere are required to avoid 
a large-angle deflection on a test particle. Hence the Debye length 
will be small in comparison to the length scale of the plasma. In 
this case, the plasma is dominated by many long-range interac¬ 
tions, rather than the short-range binary collisions of a neutral 
gas. A measure of the efficiency of the screening is the plasma 
parameter Ad- 

The plasma parameter is defined as 

Ad = (4/3)7rneA^, (Bl) 

the number of particles in a Debye sphere with radius Ad and 
centred on a single charge particle that produced the charge 
imbalance. When there are many plasma particles in a Debye 
sphere (Ad ^ 1) and long-range collective interactions are dom¬ 
inant over short-range collisions, the plasma frequency is much 
larger than the electron-ion collision frequency. Figure Bl shows 
Ad 1 for all model atmosphere structures. This indicates that 
thermal electrons interact over large distances in atmosphere of 
ultra-cool atmospheres. 

This paper has been typeset from a DT^JC file prepared 

by the author. 





Pgas [bar] 

Figure Bl. Number of particles inside of Debye sphere measure 
the efficiency of the screening given by the Debye sphere in the 
plasma. Nd ^ 1 results in the collective interactions dominate 
over short-range collisions in the gas. Top: Group 1. Middle: 
Group 2. Bottom: Group 3. 
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Table A2. First, second and third ion represent the most abundant positive ion in the high, middle and lower pressure regions respectively. 
First ion corresponds to the higher pressure (and therefore higher abundances) and visa verse. 


M/H 

Teff 


-0.6 



-0.3 



0.0 


+0.3 



1000 

K+ 

K+ 

K+ 

K+ 

K+ 

K+ 

K+ 

K+ 

K+ 

K+ 

K+ 

K+ 


Na+ 



Na+ 



Na+ 



Na+ 



1200 







K+ 

K+ 

K+ 











Na+ 






1400 







K+ 

K+ 

K+ 











Na+ 






1500 










K+ 

K+ 

K+ 











Na+ 



1600 







Na+ 

K+ 

K+ 




1800 







Na+ 

K+ 

K+ 




2000 

Na+ 

K+ 

K+ 

Na+ 

K+ 

K+ 

Na+ 

K+ 

K+ 

Na+ 

K+ 

K+ 

2200 







Na+ 

K+ 

K+ 




2400 







Na+ 

K+ 

K+ 




2600 







Mg+ 

K+ 

K+ 












Na+ 





2800 

Na+ 

K+ 

Na+ 

Na+ 

K+ 

K+ 

Na+ 

K+ 

K+ 

Na+ 

K+ 

K+ 


Ca+ 

Na+ 

K+ 

Ca+ 

Na+ 


Mg+ 

Na+ 


Mg+ 

Na+ 






K+ 









3000 







Na+ 

K+ 

Na+ 











Mg+ 

Na+ 

K+ 











H+ 








